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ABSTRACT 

Supernova measurements have become a key ingredient in current determinations of 
cosmological parameters. These sources can however be used as standard candles only 
after correcting their apparent brightness for a number of effects. In this paper we 
discuss some limitations imposed by the formalism currently used for such corrections 
and investigate the impact on cosmological constraints. We show that color corrections 
are, in general, expected to be biased. In addition, color excesses which do not add 
a significant scatter to the observed SN brightnesses affect the value of cosmological 
parameters but leave the slope of the color-luminosity relation unchanged. 

We quantify these biases in the context of the redshi ft-dependent dust ex tinction 
suggested by the recent detection of intergalactic dust by iMenard et all (I2009T) . Using 
a range of models for the opacity of the Universe as a function of redshift, we find that 
color-magnitude-stretch scaling relations are virtually insensitive to the presence of 
cosmic dust while cosmological parameters such as J7m and w are biased at the level 
of a few percent, i.e. offsets comparable to the current statistical errors. 

Future surveys will be able to limit the impact of intergalactic extinction by ob- 
serving at larger wavelengths. In addition such datasets will provide direct detections 
of intergalactic dust by cross-correlating SN colors and the density of foreground galax- 
ies, which can be used as a consistency check on the cosmic dust extinction correction. 
Alternatively, such biases could be avoided by correcting the colors of supernovae on 
an object-by-object basis with accurate photometry. 
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1 INTRODUCTION 

About a decade ago, the use of type la supernovae (SNe la) 
as standardizable candles led to the discovery of the accel- 



erating expansion of the Universe jGarnavich et al. 19981: 



Perlmutter et all Il998l ; iRiess et all ll99Sl ; ISchmidt et all 
19981 ). Since then. SNe la measurements have become a cor- 



nerstone in cosmological parameters estimation. More recent 
and upcoming surveys are aiming to push these constraints 
to higher accuracies and explore the properties of dark en- 
ergy. To this end, it is important to identify and characterize 
the potential limits and biases involved. 

The use of SNe la as standard candles requires cor- 
recting their observed magnitudes for a number of effects: a 
luminosity-light curve shape dependence (usually referred to 
as the "stretch"), a luminosity-color relation, dust extinc- 
tion and gravitational magnification. For sufficiently large 
surveys, magnification effects average out, adding only scat- 



ter to the apparent magnitudes without introducing any sig- 
nificant bias. Conversely, the effect of dust extinction is more 
critical, as it is cumulative along the line of sight. In this 
paper, we address the limits and potential biases of the pro- 
cedure currently used to correct for such effects, focusing on 
the extinction induced by intergalactic dust. 

The presence of dust on large scales around galaxies 
has been the subject of numerous papers. Theoretical stud- 
ies indicate that dust grains can be efficiently transported 
from galactic disks to the in tergalactic medium through 
wind s and radiation pressure l|Aguirrel Il999l : lAguirre et al.l 
l200ll : iBianchi fc Ferrarall2005l ). Based on estimates of the 
stellar density and metallicity as a function of redshift, 
several authors have inferred the existence of significant 
amounts of intergal actic dust and a cosmic dust density 
Oh,,,* ~ 10~ 6 -10~ 5 (lLoeb fc Haimanlll997l ; ICorasanitill2007l ; 



llnoue &: Kamavall2004 ). The amount of dust in galaxy clus- 
ters has been explored through reddening measurements 
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of background sources (e.g. Chelouche. Koester. fc BowerJ 
120071 ; iBovv. Hogg, fc Moustakasl 120081 ; iMuller et all l2008h . 
Upper limits on the the cosmic opacity were recently ob- 
tained by measuring the excess scatter seen in higher 
redshift quasar colors, dMortsell fc Goobai 
ing the Toman test l|More. Bow, fc Hogg 
combin ing constraints from lumino si ty d istances 
H(z) (|Avgoustidis. Verde, fc Jimened 120091 ) and using 
the amount of dust X-ray scattering around AGNs 



2003) 



2008) 



us- 
and 
and 



|Diikstra fc Loebll2009h. 

Recently. iMenard et al.1 (|2009l ') (M09 hereafter) cross- 
correlated the colors of about 85,000 distant quasars with 
the position of 20 million z ~ 0.3 galaxies o bserved with the 
Sloan Digital Sky Survey (|York et al.ll2000l ; SDSS). They re- 
ported a statistical detection of dust reddening up to large 
scales around galaxies. From this result they inferred the 
opacity of the Universe, calibrated at z ~ 0.3 and extrap- 
olated to higher redshifts. Their model-dependent estimate 
gives Ab(z = 1) ~ a few x 10 -2 mag. Such an opacity is 
not negligible given the precision level of current supernovae 
surveys. This Letter makes use of these results and addresses 
their impact on the estimation of cosmological parameters 
from SNe la. Our analysis does not include any effects due 
to a possible gray dust component. If present, the biases 
in cosmological parameters derived below should be seen as 
lower limits only. 



2 COSMOLOGICAL INFORMATION FROM 
TYPE IA SUPERNOVAE 

2.1 Scaling relations 

The distance modulus /i of a standard candle with intrinsic 
absolute magnitude M can be used as an estimate of the 
luminosity distance: 



ft = m - M = 5 log 10 (D L (z)/10 pc ) 



(1) 



where m is the observed magnitude. The luminosity distance 
Dl is defined by the relationship between bolometric flux S 
and luminosity L: Dl = \f L/{knS). -Dl can be related to 
the cosmology by 



(2) 



where, for a Universe containing a dark-energy component 
with density Qx and equation of state w, 



(3) 



E(z) = y / n M (l + z) 3 + fix(l + z) 3 < 1 +™) • 

The impact of dust extinction on the estimation of cosmo- 
logical parameters using the luminosity distance can be il- 
lustrated by Taylor-expanding the distance modulus from 
Equation [T] In the case of a flat Universe with a cosmologi- 
cal constant, we find 

8m ~ {—0.8, —1.3} 8Qm for a source at z = {0.5, 1} . (4) 

Therefore, the mean magnitude of standard candles must 
be controlled at a level comparable to the targeted accuracy 
of the matter density estimation. The above equations show 
that, unless corrected, the existence of a cosmic opacity of 
Ab ~ 0.05 will affect the value of Qm by a similar amount, 
translating into a ~ 15% bias for Qm = 0.3. As derived by 



IZhand (|2008h . a similar relationship can be written down for 
the shift in the equation of state parameter w: 



8m ~ 0.5 Sw 



2.2 Luminosity distance and potential biases 



(5) 



Given a sample of supernovae, constraints on cosmological 
parameters can be computed by comparing the expected 
and observed luminosity distances. To do so we define the 
X 2 -function, 



x 2 = E 



[fii - 5 log 10 (D L (zQ/10 pc)]^ 



(6) 



For each supernova, the distance modulus ^t; is calculated 
according to Equation [T] for a magnitude 



m bs,i — Snii 



(7) 



where m b s ,i is the observed magnitude (usually expressed 
in terms of the rest-frame B-band) and <5m; is the net magni- 
tude shift needed to convert the supernova from a standard- 
izable to a standard candle. Recent analyses have included 
two components to describe 8m: 

(i) a correction term to account for any color-magnitude 
relation, including intrinsic color excess, dust extinction and 
reddening, etc. 

(ii) a stretch correction to include an intrinsic dependence 
between the brightness and duration of SNe la 

The latter term is not of interest in the present analysis; it 
will be included in all calculations but not discussed explic- 
itly. 

It is generally assumed that these correction terms are 
not redshift dependent. This allows one to write the distance 
modulus as 



■ m ohSii - M + a(si - 1) - P Ci 



(8) 



where a(si — 1) refers to the stretch correction, (3 Ci is the 
color-based correction, making use of the observed color ex- 
cess Ci (usually expressed as the rest-frame E(B — V) color) 
of each supernovcQ- The coefficients a and j3 are unknown 
parameters which do not depend on redshift. The color cor- 
rection term must account for both dust extinction and any 
intrinsic color-magnitude relation with a unique coefficent 
/3 = Pa- A redshift-dependent color/extinction term will 
therefore invalidate the assumption that /3 is constant with 
z, introducing a bias into the parameter estimation. 

The color-based correction makes use of the observable 
net color excess Ci of each supernova. In general, this excess 
can be the sum of several components, 



(9) 



where the fc-terms describe the intrinsic color, dust redden- 
ing due to the host galaxy, dust along the line-of-sight and 



1 Other teams have used a different parameterization to 
apply a color-based c orrection lljha. Riess. fc Kirshnerl |2007| ; 
IWood-Vasev et al.|[2007h . Nevertheless these formalisms also rely 
on the assumption of a non-evolving color correction and similar 
conclusions apply. 
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so on. In principle, each contribution d t k can be corrected 
for using an appropriate color-dimming coefficient 



Srrii 



Ci.k 



(10) 



(if attributed to dust extinction, j3u corresponds to Rb = 
Ab/E(B — V) = Rv+1, following the previously introduced 
bands). 

Let us now consider the effects of cosmic dust extinc- 
tion, which we will denote with the subscript "d". For sim- 
plicity, we will only consider its redshift dependence, i.e. 
0i,k = /3d and a t k = Cd(Zi). The brightness of standard can- 
dles will appear modified according to 5m, = /3d Cd(zi). If 
/3d differs from the best-fit value of /3o found in Equation [5] 
a redshift-dependent magnitude bias 



8m\, 



(/3d - /3o) c d {zi) 



(11) 



is introduced. The relation between /3d and /3q depends on 
the relative contribution of cosmic dust reddening to the 
overall color scatter. Note that such a bias can be either 
positive or negative. 

In order to quantify the amplitude of this potential bias, 
we use the estimate of the opacity of the Universe given 
by M09. These authors reported reddening effects on large- 
scales around z ~ 0.3 galaxies, characterized by /3d = Ry + 
1 = 4.9 ± 2.6. Based on these results, they proposed several 
model-dependent estimates of the opacity of the Universe as 
a function of redshift. Here we will consider two cases: 

(i) a high- As model, shown in Figure [1] with the dark- 
blue curve. This estimate is calibrated at the redshift of 
their measurement and then extrapolated to higher redshift 
using an evolution model based on the observed amount of 
dust in Mgll absorbers. 

(ii) a low- As model, where the cosmic dust density for the 
high-As model is damped by a factor (l + z)~ 2 . This results 
in a significant suppression of the opacity at z > 0.8. This 
model, shown with the light-blue curve, is still consistent 
with broad observational constraints and only about a factor 
~ 2 — 3 higher than the allowed lower limits. 

Rece nt analyses of SNe la point toward a value of 
/3o — 2.5 (|Kowalski et ahlbOQct ) which is derived from min- 
imizing the residuals in the Hubble diagram. According to 
Equation [TTJ a component of cosmic dust with /3d = 4.9 will 
bias the distance modulus estimate by 



Sm hia , B (z = 0.5) ~ 0.01 mag . 



(12) 



Given the simple scaling relations derived in i )2.1l this will 
in turn bias the inferred Qm value by SQm — 0.01. This 
translates into a ~ 3% bias for Qm — 0.3. Below we quantify 
this effect more accurately, using existing data. 



2.3 Application to the Union supernova sample 

We now investigate the impact of cos mic dust on a recent su - 
pernova dataset: the "Union" sample (jKowalski et al . 2008). 
These authors have combined various supernova samples, 
compiling a "clean" dataset of 307 SNe with 0.015 < z < 
1.55 which they used to infer cosmological parameter con- 
straints. 

Using Equation [6] we sample the parameter space con- 
strained by this dataset for two cosmological models, using 



A 

CO 

< 
V 

0) 



0.10 



0.01 



~ i — i — i — r 



h — i — i — r 



n — i — i — r 



n — i — i — r~ 



T 




Tolman test (More et al. 09) 
QSO color scatter (Moertsell & Goobar 05) 
SNe la + H(z) (Avgoustidis et al. 09) 
dust around galaxies (Menard et al. 09) 
dust in Mgll absorbers (Menard et al. 07) 
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Figure 1. Top: Opacity of the Universe as a function of redshift, 
in the rest- frame B band as estimated by Menard et al. (2009). 
The light blue curve shows an extrapolation of the z ~ 0.3 mea- 
surements taking into account the evolution of dust density as a 
function of redshift. The dark blue curve applies a damping term 
to that extrapolation. Bottom: The data points show the residu- 
als of the Hubble diagram with the Union sample of supernovae. 
The blue curves show the expected contribution to the scatter in 
magnitudes due to intergalactic extinction. 



an adaptive importance sampling alg orithm called Popu la- 
tion MonteCarlo (for more details, see IWraith et aQ joojfl. 
We first consider a flat ACDM Universe using only con- 
straints from the supernovae themselves. This allows us to 
estimate the constraints on {Qm, M, q, /3}, the latter three 
being the nuisance parameters intended to absorb the ef- 
fects of the standard candle normalization. Our numerical 
values are presented in Table [1] and we show the posterior 
distributions of M, /3 and Qm in Figure [2] with the gray 
histograms. Second, we consider a dark-energy model with 
constant equation of state w. In order to obtain interesting 
constr aints, we add data fro m the CMB (WMAP5 distance 
priors; iKomatsu et a l. 2008) and BAO (distance parameter 
A measured with the SDSS; Eisenstein et al.ll2005T ). The con- 
straints on the parameters {w, Qm, Qb, h, M, a, /3} are given 
in Table [T] and we show the posterior distribution for Qm , 
Qb and w and h in Figure [3] In each case, we recover results 
similar to iKowalski et alj (fioOot ). 
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Figure 2. Posterior probabilities for the supernovae absolute 
magnitude M, the color-magnitude scaling coefficient (3 and the 
matter density fi m using the Union sample of supernovae for a 
ACDM cosmology. The different models use no dust correction 
(gray histogram), a correction for the high Ag(z) model (dark 
blue) and a low Ab(z), for which the dust density is suppressed 
at high redshifts (light blue). The dashed curve uses the high 
Ab(z) and includes the measured error on /?d- Cosmic dust does 
not affect the estimate of the nuisance parameters but does im- 
pact the constraint on !2m- 
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Figure 3. Same as Figure [2] but for a uiCDM cosmology using 
constraints from SNe Ia+CMB+BAOs. 



2.3.1 Smooth dust component 

After verifying that our method reproduces the current re- 
sults, we can now turn our attention to including the effects 
of cosmic dust extinction. First, we consider the limiting 
case of a smooth cosmic dust component. Using the two es- 
timates of the opacity of the Universe introduced above and 
(for now) ignoring measurement uncertainties, we can cor- 
rect the apparent magnitude m,i and the net color Ci of each 
supernova for the expected effect of intergalactic dust extinc- 
tion as a function of supernova redshift. In principle, such 
corrections should take into account the band used for each 
observation. In order to simplify our treatment we will how- 
ever only consider the rest-frame B-band, which corresponds 
to the observed 7?-band at the average redshift, (z) ~ 0.5, 
i.e. the filter used for most of the supernovae in this sample. 
We then use these corrected magnitudes and calculate the 
constraints on cosmological parameters as done above. 

The results of this calculation are shown in Figure [5] and 
Table Q] We first observe that the posterior distributions on 
M , a and (3 are basically unchanged. The color-magnitude 
correction used to convert supernovae into standard can- 
dles is therefore not sensitive to the presence of a redshift- 
dependent cosmic extinction. This is expected as a smooth 
component has a minimal contribution to the scatter in mag- 
nitude and/or color to which (3 is sensitive. However, the 
tilt induced in the Hubble diagram by a redshift-dependent 
cosmic extinction significantly affects the estimate of the 
matter density: &Q.M — 0.02, i.e. a ~ 7% change in am- 
plitude. This corresponds to a systematic shift comparable 
to the statistical uncertainty. By comparison, estimates of 
the absolute magnitude M and (3 change by less than 10~ 2 . 
When extending the cosmological model to a iuCDM (and 
adding constraints from CMB and BAO), we obtain the re- 
sults shown in Figure [3] As with the ACDM model, the 
best-fit values of M, a and (3 remain unchanged (see Table 
[TJ, while the cosmological parameters are all biased to some 
extent. As shown in the Table, the offsets range from 0.25 
to 0.45cr. 

For a more realistic description of the properties of the 
intergalactic dust we can include the measurement error on 
/3d- To do so, we add a Gaussian prior for /3d to the \ 2 with 
mean 4.9 and rms 2.6 as indicated by M09. The results for 
the high- Ab(z) model are shown in Figures [2] & [3] using 
the dark-blue dashed lines. The uncertainty in (3d propa- 
gates into the best-fit parameter uncertainties, yielding an 
increase in the width and skewness of the posterior distri- 
butions. 

These results illustrate that with the current formalism 
being used to measure the luminosity distance from SNe la, 
the color-magnitude correction term does not properly ac- 
count for a redshift-dependent cosmic dust extinction. The 
slope /3 of the color-luminosity relation is not sensitive to 
such a reddening component. As a result, the formalism can 
only correct for a fraction of the extinction effect, as in- 
dicated in Eq. 1111 Thus, instead of being absorbed by the 
nuisance parameters, such a component of dust impacts the 



2 In Equation [6] the error er^ for each supernova contains the co- 
variance of rdoba i,Si, Ci (which depends on a and /3), an intrinsic 
absolute magnitude scatter of 0.15 and an additional uncertainty 
from peculiar velocities of 300 km s — 1 . 
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Parameter 


No Correction 


High A B 
Pd = 4-9 


High A B 
/3 d = 4.9 ± 2.6 


Low Ab 
/3 d = 4.9 






ACDM: U M 


2gl +o.032 
u '^ a -0.030 


0.308lg;^ (0.55o-) 


0-308«;[J39 ( .55ct) 


0.304«;033 (o.42a) 


Parameter 


all models 


wCDM: n h 


n 0457+ - 002 

u - u ^ ' -0.002 


0.046±°;°™ (0.35a) 


0.04B±g;gg (0.25ct) 


0-045«;^ ( .25ct) 


M 

a 


-19.31 ±0.03 
1.37 ±0.13 


h 


U.DMO_ Q17 


0.687±S".oi? (0-45o-) 


0-688t°;°?° (0.40ct) 


0.688to gi7 (0.40cr) 


13 


2.45 ± 0.12 




n 97 o+0.017 
u - z '°-0.016 


0-279±g;gi2 (0.36a) 


0.278to.ol? (0.30cr) 


o-^toTie (°- 30ff ) 






— w 


0.968ir o 6 6 ! 


0-940± -°tl (0.43a) 


0.944i°;°g (0.37a) 


0.944±°;°«| (0.37a) 







Table 1. Constraints on the fitting parameters of Equation [8] for both a ACDM cosmology (using SNe la) and a toCDM cosmology 
(using SNe Ia+CMB+BAO). Three estimates of the opacity of the Universe are considered (see text). Note that in all cases, we obtain 
the same best fit values for a, (3, M. In other words, the current parameterization is not sensitive to a component of cosmic dust. We 
find that such a component can bias the estimate of £Im and w by 0.3 to 0.5 <r with current data. 



value of cosmological parameters, resulting in offsets com- 
parable to the current statistical errors on Qm and w. 

2.3.2 Clumpy dust component 

Since dust is expected to originate from galaxies, the ex- 
tinction effects due to intergalactic dust should be related 
to the underlying density field. Hence, the dust contribution 
to the scatter in supernovae colors/magnitudes should re- 
flect this correlation. Using a toy model we can investigate 
the expected impact of this effect on cosmological parameter 
estimation. 

The number of halos intercepted by an average line-of- 
sight is given by 

N S = f\n^±*f Az > (13) 

where a and n are the cross-section and number density 
of the galaxies responsible for most of the dust extinction. 
Assuming that the amount of dust is roughly proportional 
to luminosity and metallicity, M09 claim that most of the 
dust originates from galaxies with L ~ L* /2. They use n ~ 
0.037 /i 3 Mpc" 3 and a = vrr^, with r v ~ 110/i _1 kpc. 

For a Poisson distribution, the level of fluctuation is 
simply ~ 1 / *J N g . An upper limit on the extinction scatter 
due to intergalactic dust is then given by Ab(^)/v N g (z). 
This quantity is shown in the lower panel of Figure [1] for the 
two models of dust considered above. At z = 1, the average 
number of halos intercepted by a line-of-sight is of order ten 
and the scatter in Ab is about 0.03 mag. 

Comparing this quantity to the scatter in observed su- 
pernova magnitudes from the Union sample (data points), 
we see that the the level of scatter due to cosmic extinction 
ends up being significantly lower than the observed scatter in 
SN magnitudes. Therefore, this contribution is not expected 
to substantially affect the above constraints on cosmological 
parameters and leaves our previous conclusions unchanged. 



3 DISCUSSION AND PROSPECTS 

The current formalism used to correct the brightness of su- 
pernovae for color-magnitude trends assumes that all color- 



magnitude corrections can be parametrized by a unique 
(and therefore redshift independent) coefficient, /3. We have 
shown that in the presence of several sources of reddening 
this assumption may lead to a biased estimate of the intrin- 
sic brightness of SNe and therefore their distance modulus. 
In addition, color excesses which do not add a significant 
scatter to the observed SNe la brightnesses end up affecting 
the value of cosmological parameters but leave the best-fit 
color-scaling parameter /3 unchanged. 

We have considered the effects of the intergalactic dust 
observed by M09 and applied a range of extinction correc- 
tions to the Union supernovae sample. We have shown that 
the reddening-based correction used in recent works is not 
sensitive to such an extinction contribution. As a result, it 
only corrects for a fraction of the extinction effect, and biases 
estimates of cosmological parameters. Considering several 
scenarios of extinction corrections as a function of redshift 
we found biases in SIm and mat a level ranging from 0.3 to 
0.5a of current statistical errors. 

While our analysis indicates a systematic limitation in 
the current formalism for standardizing SNe la, we remind 
the reader that the parameter estimates presented above 
rely on the value of the opacity of the Universe which re- 
mains poorly constrained, especially at high redshift. The 
quantitative results are based on reddening measurements 
around galaxies at z ~ 0.3, then extrapolated. Our analysis 
is therefore not meant to provide more robust cosmological 
parameter estimations, but only illustrates the fact that the 
statistical errors achieved by recent surveys are becoming 
comparable to systematic biasing due to extinction. Below 
we suggest a few means by which future surveys can reduce 
their exposure to this particular bias. 



3.1 Minimizing the effect of dust extinction 

Dust extinction effects weaken significantly when consider- 
ing larger wavelengths. One can therefore minimize the am- 
plitude of the bias discussed above by observing and estimat- 
ing supernova magnitudes in redder bands. Planned surveys 
will satisfy this criterion: 

(i) JDEM will provide us with observations of supernovae 
in rest-frame I band. For dust with Rv = 3, the extinction 
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effects are expected to be about two times weaker than those 
in B band. According to Equation the bias in w will de- 
crease to about 1%. More detailed calculations applying the 
extinction corrections to simulated JDEM data give about 
0.5% (D. Rubin, private communica tion) . 

(ii) For LSST (jlvezic et al.l I2008D the reddest filters (z 
and Y) should be sufficient for k-correcting to rest-frame R 
band, which should yield roughly a 1.5% systematic error 
on w from intergalactic dust using the current formalism. 

In addition, if the photometry is accurate enough, it 
becomes possible to correct for reddening effects object-by- 
object. With photometry in three or more passbands and a 
supernova template, it is possible to estimate Rv and E(B — 
V) for each object. Such an approach no longer requires 
the assumption of a unique slope /3 of the supernova color- 
luminosity r elation and might be favo red in the future. As 
pointed out IWood-Vasev et ail l|2007l ). there are currently 
only a few SNe la in the literature with the requisite high- 
precision photometry extending from the rest-frame UV to 
the near-IR, at z > 0.2. 



E B -v 10~ 2 and N g ~ 7, giving N s ~ 130. Such detec- 
tions might therefore be within reach, given the expected 
yields of the next generation of supernova surveys. Hence, 
the contribution of cosmic dust extinction can be corrected 
for, on average, as done in the present analysis. Note that 
Ns oc <t(c) 2 . Therefore, better photometric accuracy would 
significantly reduce the number of required supernovae to 
detect the effects of intergalactic reddening. 

We are also pointing out that once the mean redden- 
ing per galaxy has been estimated, it is possible to correct 
the supernova magnitudes object-by-object, by subtracting 
the expected amount of extinction given the number of fore- 
ground galaxies in the vicinity of each supernova. If the am- 
plitude of intergalactic reddening is properly estimated, the 
correlation between foreground galaxy density and corrected 
supernova magnitudes should vanish. This provides a con- 
sistency check for the correction of cosmic dust. Such a test 
emphasizes the necessary interdependence of the large-scale 
structure and synoptic aspects of the next-generation sur- 
veys. 



3.2 Detecting and correcting for intergalactic 
dust extinction with SNe la 

Similarly to search es for magnification effects 
jMenard fc Dalall l2005h . the contribution of dust ex- 
tinction associated with foreground galaxies can in principle 
be detected using the supernovae themselves, by cross- 
correlating their observed colors with the density field of 
foreground galaxies: 

(Sc) g (e) = ([c(^-{c}]s s (<p + e)), (14) 

where 

8 g (9) = N g (6)/(N g ) - 1 (15) 

and iVg is the density of foreground galaxies measured within 
some aperture of size 8 around a given supernova. Note that 
including higher-redshift galaxies would not bias the esti- 
mate but only add noise. Compared to the results presented 
in this paper, the advantage of such a measurement is to 
directly include the full redshift range from which dust ex- 
tinction effects arise. Similarly, measuring this effect as a 
function of wavelength provides constraints on the coeffi- 
cient /3d relevant to the considered supernovae sample. 

We now estimate the number of supernovae required 
for such a measurement. The above correlation is not sen- 
sitive to the mean reddening along lines-of-sight but only 
to fluctuations. The typical reddening excess due to inter- 
galactic dust is given by Eb-v(z) a(N g ), where N g is de- 
fined in Equation [13] Such a change will be induced to the 
mean color of supernovae lying behind the corresponding 
over/under-densities. The error on the estimate of the mean 
color of A^s supernovae is given by ai c \ = <r(c) / V Ns, where 
cr(c) is the scatter of observed supernova colors. Detecting 
such a quantity at the v er-level requires 

^<Bb-yWVW. (16) 
V Ns 

Atz = 0.5, (Eb-v) ^5x 10~ 3 and N g ~ 3. For <j(c) ~ 0.1 
(see Figure [TJ, we find that Ns ~ 1200 supernovae are 
required to detect this effect at the 3a level. At z = 1, 
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